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AMENDMENTS TO THE CLAIMS : XIF-RSC'Cl PGT/PTO .0 2 MAY 20D5 

This listing of claims will replace all prior 
versions, and listings of claims in the application: 

LISTING OF CLAIMS : 

1 . (original) A digital simulation method of a non-linear interaction between an 
excitation source and a wave in a resonator, by means of digital signal calculation tools 
based on equations the solution of which corresponds to the physical event of a 
phenomenon to be simulated which can translated, at each time and at each point of the 
resonator, by a linear relation between two variables representative of the effect and of the 
cause of said phenomenon to be simulated, a method wherein the impedance or 
admittance equation is transcribed directly into a digital model form enabling to realise a 
non-linear interaction between the two variables of the impedance or admittance relation. 

2. (original) A method according to claim 1, for simulating an 
oscillating phenomenon, wherein the digital model comprises, on the one 
hand at least one linear part (3), representing the input impedance or 
admittance of the resonator, and, on the other hand, one non-linear part 
(2) modelization the role of the excitation source (22) of the phenomenon 
to be simulated. 

3. (currently amended) A simulation method according to c l aims 1 or 2 claim 1 . for 
real time digital synthesis of an oscillating phenomenon, wherein, from a system of 
equations between at least two variables representative of the behaviour of a wave in the 
resonator, an expression of the input impedance or admittance of the resonator is 
established in the form of a linear filter including delays, without any decomposition into 
two-way waves, in order to realise at least one linear part (3) of the model. 

4. (original) A method according to claim 3, wherein the linear part (3) of the model 
is coupled with one non-linear part (2) involving the evolution of the non-linearity as 
expressed between the two variables of the input impedance or admittance relation of the 
resonator. 

5. (original) A method according to claim 4, wherein the linear part (3) of the digital 
simulation model of the impedance or admittance equation in based on to two elementary 
waveguides fulfilling a transfer function between the two variables of the impedance or 
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admittance relation. 

6. (original) A method according to claim 5, wherein the linear part (3) with two 
waveguides of the model is coupled with a loop connecting the output to the input of said 
linear part (3) and comprising a function (21) involving the non-linearity as expressed 
physically. 

7. (original) A method according to claim 6, wherein the model is driven by at least 
two parameters representative of the non-linear physical interaction between the source 
and the resonator, by means of a loop connecting the output to the input of the linear part 
(3) and comprising a non-linear function (21) playing the part of an excitation source for the 
resonator. 

8. (currently amended) A method according to one of th e proviouc claims claim 1 . 
for real-time sound synthesis of a musical instrument comprising, at least, one excitation 
source with non-linear characteristics and a linear resonator, the sound produced by the 
instrument resulting from a coupling, between the excitation source and the resonator, 
expressed at least by a linear impedance or admittance relation and a non-linear relation 
between two physical variables representative of the effect and of the cause of the sound 
produced, a method where the sound produced by the instrument is simulated, in real 
time, by modelization the physical phenomena governing the operation of the instrument, 
wherein, to realise said physical modelization, the impedance or admittance linear relation 
is expressed directly and digitally between two physical variables representative of the 
cause and of the effect of the phenomenon to be simulated and said impedance or 
admittance relation in digital fonn is associated with the non-linear relation between the 
same variables. 

9. (original) A method according to claim 8 for synthesis of the sound of an 
instrument with complex resonator, wherein the resonator is decomposed into a series of 
successive elements and wherein the impedance or admittance relations corresponding 
respectively to each element of the resonator are calculated and combined in order to 
obtain a global impedance corresponding to the geometry of the resonator. 

1 0. (currently amended) A method according to claim 8 [[or claim 9]], wherein, for 
real-time synthesis of the sound produced by a wind instrument, the two variables of the 
impedance relation are the acoustic pressure (pe) and flow (Ue) at the input of the 
resonator. 

11. (original) A method according to claim 10, wherein, for an open-ended 
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cylindrical resonator, the linear part (3) of the digital transcription model of the impedance 
equation is the sum of two elementary waveguides having as excitation source the flow 
(Ue) at the input of the resonator, and fulfils the transfer function : 



' Ue(cD) 1 + exp(-2ik((o)L) 1 + exp(-2ik(<o)L) 
wherein: 

CO is the angular frequency of the wave ; 
Ze(a>) is the input impedance of the resonator, 

Pe(a)) and Ue(a)) are the Fourier transforms of the dimensionless 
values of the pressure and of the flow at the input of the resonator ; 

k(o>) is a function of the angular frequency which depends on the 
phenomenon to be simulated ; 

L is the length of the resonator, 

1 2. (original) A method according to claim 1 1 , wherein each of the two waveguides 
involves a filter having as a transfer function: 

- F((o)' = - exp(-2ik((D)L) 

and representing a two-way travel of a wave, with a sign change at the open end of 
the resonator, each waveguide corresponding to a term of the Impedance equation. 

13. (original) A method according to claim 12, wherein the model is driven by the 
length (L) of the resonator and at least two parameters (i;, y) representative of the non- 
linear physical interaction between the pressure (pe) and the flow (u©) at the input of the 
resonator, by means of a loop connecting the output to the input of the linear part (3) and 
comprising a non-linear function (21 ) as an excitation source for the resonator. 

14. (original) A method according to claim 1 1 , wherein the non-linear function has 
the pressure and the displacement of the vibration formation member as input parameters, 
and is controlled by at least two parameters simulating a player playing. 

15. (original) A method according to claim 14, wherein the playing 
parameters for controlling the non-linear function are: 

a parameter q characteristic of the mouthpiece and of the action of 
the player on the vibration formation member. 



Ze(co) = 



Pe((o) 



exp(-2lk(co)L) 



(12) 
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a parameter y representative of the pressure applied to the vibration 
formation member. 

16. (currently amended) A method acxx)rding to on e of th e claims 1 0 to 1 5 claim 
10 . wherein, for real-time synthesis of the sound to be simulated, a formulation is realised 
in the time domain of the angular frequency response of the impedance of the resonator, 
by approximation of the losses represented by the filter by means of an approximated 
digital filter. 

17. (original) A method according to claim 16, wherein, to express the angular 
frequency response of the impedance of the resonator, a one pole digital filter is used of 
the form: 



^ b, exp(-2i^D) 

1 - a^ exp(-iTn) 

(13) 

wherein : 

, /e being the sampling frequency, 

c is the pure delay corresponding to an away or return travel 
of the wave in the resonator, 

- the coefficients bo and a1 are expressed in relation to the physical 

parameters SO that F((o)^ = r(txj) for a value ^ _ . , 

' ' • I (o1 of the angular frequency 

corresponding to the fundamental playing frequency and another value co2 

corresponding to a harmonic, 

and the following differential equation is derived: 

Pe(n) = UeCn) - a,Ue(n - 1) - boU,(n - 2D) + a,p,(n - 1) - boPe(n - 2D) (16) 



18. (original) A method according to claim 17, wherein the 
coefficients bo and a1 are obtained by solving the equation system: 
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|F((Oi)^f (1 + - 2ai cos(tni)) = bo 



|F(«)2)^r(1 + - 2ai cosCtUa)) = 



with |F(a))^|^ = exp(-2ac^L), said coefficients being given by the 



formulae: 



a, = 



A1 - A2 - V(A, - A,)^ - (F, 



^j2F,F,{c, - c^KA, - A, - /(A, - A,f - (F, - F,f) 



° F - F 



2 



wherein 

= cos(wJ,C2 = cos(tn2),Fi = |F((Oi)^| .Fg = |F(g)2)^| , 
= F^C-|, A2 — ^2^2 

1 9. (original) A method according to claim 1 8, for simulating a cylindrical resonator 
instrument, from a physical modelizatlon governed by the system of equations : 



1 <^^-W,3L^^^^t)^^p^^t) 



0)2 dt2 ^r dt 



(with the sign + for a reed and the sign - for the lips) 



Pe(a)) = itan — -i acco^'^L \J^{(o) 

I c 2 J 



Ue(t) = ^(1 - sign(Y - x(t) - 1))sign(Y ^ Pe(t))C(1 - Y + x(t))J|Y - p,(t)| 

^ wherein 
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or is the resonance frequency and qr is the quality factor of the reed or of 
the lips, wherein that said system of equations is solved in the time 
domain from an equivalent sampled formulation of the angular frequency 
response of the displacement of the reed or of the lips and of the 
impedance relation which is translated by the system of equations : 

x(n)=biaPe(n-1 )+aiaX(n-1 )+a2aX(n-2) 

(18) 

Pe(n)=Ue(n)-aiUo (n-1 )-boUe(n-2D)+aiPe(n-1 )boPe(n-2D) (19) 

Ue(n) = ^ (1 - sign(Y - x(n) - 1))sign(y - Pe(n))^n)- y + x(n))^|Y - p,(n)l (20) 

said equations being used sequentially by grouping the terms not depending on the 
time sample n, in order to calculate in succession : 

x(n)=biaPe(n-1 )+aiaX(n-1 )+a2aX(n-2) 

(21) 

V=-aiUe(n-1 )-boUe(n-2D)+aiPe(n-1 )-boPe(n-2D) 

(22) 



W = |(1 - sign(y - x(n) - 1))Q) - y + x(n)) 

(23) 

1 . 



Ue (n) = ^ sign(Y - V)(-bCo + w7(bCoW)^ + 4|y - V| (24) 



Pe(n)=boCoUe(n)+V 



(25) 



20. (original) A method according to claim 19, for more realistic simulation of the 
produced sound, wherein, by neglecting the radiation, the external pressure is expressed 
as the time derivation of the outgoing flow, in the fomri: 
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Pext(t) = |(Pe(t) + Ujt)) 

(26) 

and is calculated, at each sampled time (n), by differential between the sums of the 
intemal pressure pe and of the flow ue. respectively at the time (n) and at time (n-1 ). 

21. (original) A method according to claim 19 for simulating a multimode reed 
instrument, wherein the calculation of the acoustic pressure and of the flow at the 
mouthpiece is performed by sequential resolution of a system of equations wherein the 
displacement of the reed at each time(n) is in the form: 

x(n)=baiPe(n-1 )+ba2Pe(n-2)+baDiPe(n-Da-1 ) +aaix(n-1 ) +aa2x(n-2)+aaDx(n- 
Da)+aaDix(n-Da-1 ) 

the coefficients aal, aa2, aaD2, aaD1 being defined by: 



^ .fe(1±MZ_P ao-^^iP^ an=-b a..- ^a(fe -P) 
aal f aa2 I oaD — Da, aaD1 ' 

and the coefficients bal, ba2, baDI by : 




(28) 



by setting : = . , A2=-cD,q, and A3=A2Aiqr, 

the following equations being the same as for a single mode reed. 

22. (currently amended) A method according to one of tho c l a i ms 1 9 
to 21 claim 19 , wherein that a model is prepared for a cylindrical 
resonator with terminal impedance from the basic model corresponding to 
a cylindrical resonator and being the sum of two waveguides involving 
each a filter having as a transfer function -F(co)^=-exp(-2ik(co)L), while 
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replacing the expression exp(-2ik(cD)L) by the expression 
R((D)exp(-2ik(co)L). wherein R(a)) = 

Zc +Zs(^o) 
Zc being the characteristic impedance 

-^etZ^ the output impedance 

23. (currently amended) A method according to on e of th e cla i ms 19 to 22 claim 
19 . wherein, from the impedance model for cylindrical resonator instrument and the 
associated differential equations, other more complex impedance models are built for 
simulating oscillating phenomena produced by a resonator of any shape by combining 
impedance elements in parallel or in series and by using digital approximations for an 
explicit use of the physical variables involved in the production of said oscillating 
phenomena and a more flexible control of the result of the simulation. 

24. (cun^ently amended) A method according to on e of the c l a i ms 19 to 23 claim 
23. wherein that, from the basic model of a cylindrical resonator wherein the angular 
frequency response of the displacement of the reed or of the lips is translated by a system 
of differential equations providing, at each time(n), the displacement x(n) the pressure 
pe(n) and the flow ue(n) at the input of the resonator, a model for a conical resonator is 
built wherein the equation of the pressure is in the form: 

Pe(n) = bCoUe(n) +bciUe(n-1)+bC2Ue(n-2)+bCDUe(n-2D)+bCDiUe(n-2D-1) 

+ aciPe(n-1 )+ac2Pe(n-2)+acDPe(n-2D)+acDiPe(n-2D-1 ) (33) 

wherein the coefficients bcO, bcl, bc2, bcD and bcDI are defined by 



1 3- + 1 . a. . . bo 

bCo= — . bci=- — , bC2=— ^, bCD=- 7^^. bCDi=-r- 
Gp Gp Gp Gp Gp 

and the coefficients aci, ac2, aco and acoi are defined by : 



ac^ = - 



aiGp 



acj = 



G„ 



acp = 



boG^ 



acp^ = bo 
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1 1 
by noting : Gp = 1 + and Gm = 1 



2f e ^ 2f 3 ^ 

c c 

25. (currently amended) A method according to on e of th e claims 1 9 to 23 claim 24 

wherein that, from the basic model for a cylindrical resonator, a model for a short resonator 

having a length I is built, by an approximation of the impedance according to the 

expression: 

Z,((d)= i tan (k((D)l) = G(co) + icoH(a)) (34) 



1 _ expf- ac^ I 

wherein G(co) = f j=4 and H ((o)=-(l - G((d)). 

1 + exp 



- acJ- I 



26. (currently amended) A method according to cla i ms 23, 2A, ciaim. 
25 i n comb i nat i on, for simulating a wind instrument, wherein the 
mouthpiece or the bill is modelled by a Helmholtz resonator comprising a 
hemispheric cavity coupled with a short cylindrical pipe and a main 
resonator with a conical pipe, the input impedance of the resonator 
assembly which may be expressed as: 

J_ 



Ze(C0) = 



. V 1 

ICO- 



P^' ico^ 



iZ,(ki(co)Li) + Z2 



1 + 



ICO — 
C 



itan(k2 ((0)1^2) 



wherein V=— tiR^ is the volume of the hemispheric cavity, L1 is the 
6 

length of the short pipe, L2 is the length of the conical pipe, Z1 and Z2 
are the characteristic impedances of both pipes which depend on their 
radii, k1(cD) and k2((o) take into account the losses and the radius R1 and 
R2 of each pipe, and that, from the basic model for cylindrical resonator, 
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from its extensions to the conical pipe and to the short pipe, a resonator 
model is prepared by expressing the pressure at the mouthpiece or in the 
bill by the differential equation : 

Pe(n)='£bc,U3(n - k) + ZbCD,U3(n - k - 2D) (36) 

k=0 k-0 
k=4 k=3 

+ XaCkPe(n - k) + X^^DkPeCn - k - 2D) 

k=1 k«0 

27. (currently amended) A method according to on e of th e c l aims 1 to 9 claim 1 , for 
the simulation of an oscillating phenomenon wherein both physical variables of the linear 
relation are the strength applied to one point of a mechanical system such as a string 
generating vibrations and the speed at this point, wherein that the admittance is 
expressed, at this point, in the form of a combination of the admittances of each part of the 
string, on both sides of said point, each mechanical admittance being obtained from the 
basic model describing the acoustic impedance of a cylindrical pipe resonator, by 
expressing the speed at the point considered of the string in relation to the strength applied 
to this point, where the filter F(co) of the basic model may be expressed from a bending 
wave propagation model in a having a stiffeness. 

28. (currently amended) A digital device to implement the method according to ©ne 
of tho claimc 8 to 27 claim 8 , for the simulation of a musical instrument generating a sound 
resulting from a coupling between a linear resonator and an excitation source with a non- 
linear characteristic, which can be expressed at least by a linear impedance or admittance 
relation and a non-linear relation between two variables representative of the effect and of 
the cause of the produced sound, said simulation device comprising a control element I 
including at least one gestural sensor 1 transforming the actions of a player into control 
parameters, a modelization element II including one non-linear part (2) associated with a 
linear part (3) and an element creating the synthesised sound III, wherein that the linear 
part (3) includes a computing block (31 ) driven by the length (L) of the resonator, having as 
an input parameter a signal representative of one of the variables, cause or effect, 
computed by the non-linear part (2), and the transfer function of which is the input 
impedance or admittance of the resonator, in that the non-linear part (2) implements a non- 
linear function (21) driven by at least two control parameters and having as input 
parameters a signal representative of the other variable, cause or effect, computed by the 
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linear part (3) and a signal modelization the role of the excitation source, the linear part (3) 
being thus coupled with the non-linear part (2) in a closed loop, and in that the element 
creating the sound III computes a sound signal from signals representative of the cause 
and of the effect of the sound to be simulated, emitted respectively by the linear part (3) 
and the non-linear part (2). 

29. (new) A method according to claim 1 9, wherein that, from the basic model of a 
cylindrical resonator wherein the angular frequency response of the displacement of the 
reed or of the lips is translated by a system of differential equations providing, at each 
time(n), the displacement x(n) the pressure pe(n) and the flow ue(n) at the input of the 
resonator, a model for a conical resonator is built wherein the equation of the pressure is in 
the form: 

Pe(n) = bCoUe(n) +bCiUe(n-1)+bC2Ue(n-2)+bCDUe(n-2D)+bCDiUe(n-2D-1) 

+ aciPe(n-1 )+ac2Pe(n-2)+acDPe(n-2D)+acDiPe(n-2D-1 ) (33) 

wherein the coefficients bcO, bcl , bc2, bcD and bcD1 are defined by 



bco=— . bci=--^:^ — , bC2=-^. bCo=-~, bCDi=-^ 
Gp Gp Gp °P 

and the coefficients aci, ac2, aco and acoi are defined by : 



ac- = 



acg = 



aCn = - 



boG, 



aCoi = bo 



by noting 



Gp=1 + 



1 



2fe^ 
^ c 



and Gm = 1- 



2f. 



30. (new) A method according to claim 1 9 wherein that, from the basic model for a 
cylindrical resonator, a model for a short resonator having a length I is built, by an 
approximation of the impedance according to the expression: 
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Z|(co)= i tan (k((o)l) = G(<o) + ia)H(a)) 



(34) 



wherein G(co) = 



- acJ— I 
V2 



1 - exp - acJ— I 



1 + exp 



andH(a))=-(l-G(co)). 
c 
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